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Abstract. The dynamical mean-field theory (DMFT) is a widely applicable approx- 
imation scheme for the investigation of correlated quantum many-particle systems 
on a lattice, e.g., electrons in solids and cold atoms in optical lattices. In particu- 
lar, the combination of the DMFT with conventional methods for the calculation of 
electronic band structures has led to a powerful numerical approach which allows 
one to explore the properties of correlated materials. In this introductory article we 
discuss the foundations of the DMFT, derive the underlying self-consistency equa- 
tions, and present several applications which have provided important insights into 
the properties of correlated matter. 



1.1 Motivation 

1.1.1 Electronic Correlations 

Already in 1937, at the outset of modern solid state physics, de Boer and 
Verwey [1] drew attention to the surprising properties of materials with in- 
completely filled 3d-bands. This observation prompted Mott and Peierls [2] 
to discuss the interaction between the electrons. Ever since transition metal 
oxides (TMOs) were investigated intensively [3]. It is now well-known that 
in many materials with partially filled electron shells, such as the 3d transi- 
tion metals V and Ni and their oxides, or 4/ rare-earth metals such as Ce, 
electrons occupy narrow orbitals. The spatial confinement enhances the effect 
of the Coulomb interaction between the electrons, making them "strongly 
correlated". Correlation effects can lead to profound quantitative and qual- 
itative changes of the physical properties of electronic systems as compared 
to non-interacting particles. In particular, they often respond very strongly 
to changes in external parameters. This is expressed by large renormaliza- 
tions of the response functions of the system, e.g., of the spin susceptibility 
and the charge compressibility. In particular, the interplay between the spin, 
charge and orbital degrees of freedom of the correlated d and / electrons and 
with the lattice degrees of freedom leads to an amazing multitude of order- 
ing phenomena and other fascinating properties, including high temperature 
superconductivity, colossal magnetoresistance and Mott metal-insulator tran- 
sitions [3]. 
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1.1.2 The Hubbard Model 

The simplest microscopic model describing interacting electrons in a solid is 
the one- band, spin- 1/2 Hubbard model [4-6] where the interaction between 
the electrons is assumed to be so strongly screened that it is purely local. 
More generally the Hubbard model applies to lattice fermions with a point 
interaction, such as ultra-cold fermionic atoms in optical lattices where the 
interaction is indeed extremely short ranged. The Hamiltonian consists of two 
terms, the kinetic energy Hq and the interaction energy Hi (here and in the 
following operators are denoted by a hat): 

H = H + Hi (1.1a) 

i,j cr k.a 

Hi = U^n i ^h ii , (1.1c) 

i 

where cf a (c i(7 ) are creation (annihilation) operators of fermions with spin a at 
site Ri (for simplicity denoted by i), and h ia — cf a c ia . The Fourier transform 
of the kinetic energy in (1.1b), where t,j is the amplitude for hopping between 
sites i and j, involves the dispersion e/, and the momentum distribution op- 
erator n ka . This model provides the basis for most of the theoretical research 
on correlated electrons during the last decades. 

The Hubbard model describes an interacting many-body system which 
cannot be solved analytically, except in dimension d = 1 for nearest-neighbor 
hopping [7]. In the absence of exact solutions there is clearly a great need 
for reliable, controlled approximation schemes for this model. However, such 
approximations are not easy to construct as the following observation shows. 
When viewed as a function of time a given site of the lattice will sometimes 
be empty, singly occupied or doubly occupied. For strong repulsion U dou- 
ble occupations are energetically very unfavorable and are therefore strongly 
suppressed, implying (n^n^) 7^ (n^)(n^). Therefore approximation schemes 
based on the factorization of the interaction term, e.g., Hartree-Fock-type 
mean-field theories, are generally insufficient to explain the physics of electrons 
in their paramagnetic phase beyond the limit of weak interactions. This is due 
to the fact that in such approximations the interaction is described only as 
an average, static potential, whereby correlations, i.e., dynamical many-body 
effects due to the interaction of individual electrons, are excluded from the 
beginning. Hence correlation phenomena such as the Mott-Hubbard metal- 
insulator transition cannot be described by such approximations. This clearly 
shows the need for comprehensive approximation schemes, which are applica- 
ble for all values of the input parameters, e.g., coupling parameters and tem- 
perature, diagrammatically controlled, and thermodynamically consistent [8]. 



1.1.3 Construction of Comprehensive Mean-Field Theories for 
Many- Particle Models 

There exists a well-established branch of approximation techniques which 
makes use of the simplifications that occur when some parameter is taken to 
be large (in fact, infinite), e.g., the length of the spins S, the spin degeneracy 
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N, the spatial dimension d, or the coordination number Z, i.e., the number of 
nearest neighbors of a lattice site. 3 Investigations in this limit, supplemented 
if possible by an expansion in the inverse of the large parameter, often provide 
valuable insights into the fundamental properties of a system even when the 
inverse parameter is not very small. 

One of the best-known mean-field theories in many-body physics is the 
Weiss molecular-field theory for the Ising model [10]. It is a prototypical single- 
site mean-field theory which becomes exact for infinite-range interaction, as 
well as in the limit of the coordination number Z — > oo or the dimension 
d — > oo. In the latter case \jZ or 1/d is a small parameter which can be 
used to improve the mean-field theory systematically. This mean-field theory 
is comprehensive in the sense discussed above. Namely, it contains no unphys- 
ical singularities, is applicable for all values of the input parameters, i.e., the 
coupling parameter, magnetic field, and temperature, and is diagrammatically 
controlled [11]. 

Itinerant quantum mechanical models such as the Hubbard model and its 
generalizations are much more complicated than classical, Ising-type models. 
Generally there do not even exist semiclassical approximations for such mod- 
els that might serve as a starting point for further investigations. Under such 
circumstances the construction of a mean-field theory with the comprehensive 
properties of the Weiss molecular field theory for the Ising model will necessar- 
ily be much more complicated, too. Here the limit of high spatial dimensions 
d or coordination number Z has again been extremely useful since it provides 
the basis for the construction of a comprehensive dynamical mean-field theory 
(DMFT) for lattice fermions. 

In this article we will first discuss (sec. 1.2) the limit of high spatial di- 
mensions d for lattice fermions, the scaling of the hopping amplitude which is 
necessary to obtain a meaningful limit d — > oo as well as the simplifications 
of the many-body perturbation theory occurring in this limit. In sec. 1.3 the 
self-consistency equations obtained in the limit d — > oo are derived which pro- 
vide the basis for the DMFT for correlated lattice fermions. An example for 
the many insights gained by the DMFT is the Mott-Hubbard metal-insulator 
transition discussed in sec. 1.4. The application of the DMFT to real corre- 
lated materials is described in sec. 1.5. Brief introductions to the DMFT for 
correlated systems in the presence of disorder (sec. 1.6), correlated bosons in 
optical lattices (sec. 1.7), and systems in non-equilibrium (sec. 1.8) are also 
given. In sec. 1.9 a summary and outlook is presented. 



3 The coordination number Z is determined by the dimension d and the lattice 
structure. Already in d — 3 the coordination number can be quite large, e.g., 
Z = 6 for a simple cubic lattice, Z — 8 for a bec lattice and Z = 12 for an 
fcc-lattice, making its inverse, 1/Z, rather small. It is then natural to consider 
the limit Z — > oo to simplify the problem. For a hypercubic lattice, obtained 
by generalizing the simple cubic lattice in d = 3 to arbitrary dimensions, one 
has Z — 2d. The limit d — > oo is then equivalent to Z — > oo. Several standard 
approximation schemes which are commonly used to explain experimental results 
in dimension d = 3 are exact only in d, Z — oo [9] . 
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1.2 Lattice Fermions in the Limit of High Dimensions 
1.2.1 Scaling of the Hopping Amplitude 

We consider the kinetic energy term (1.1a) since the interaction term is purely 
local and is thereby completely independent of the lattice structure and the 
dimension. For hopping between nearest-neighbor (NN) sites i and j with 
amplitude tij = — t on a <i-dimcnsional hypercubic lattice with lattice spacing 
a, the dispersion e k is given by 

d 

e k = -2tJ2 cos(k n a). (1.2) 

ra=l 

The density of states (DOS) corresponding to e k is 

N d (Lu)=Y,^-e k ), (1.3) 
k 

which is the probability density for finding 4 ui = e k for a random choice of 
k = (ki, . . . ,kd). If the ki are chosen randomly e k in (1.2) is the sum of 
(independent) random numbers — 2t cos ki. The central limit theorem then 
implies [12] that in the limit d — »■ oo the DOS is given by a Gaussian 
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(1.4) 



Unless t is scaled properly with d this DOS will become arbitrarily broad and 
featureless for d — > oo. Clearly only the scaling 

t — > -^=, t* = const., (1.5) 
V d 

("quantum scaling") yields a non-trivial limit d — >• oo for the DOS [13, 12]. 

The interaction term in (1.1) is seen to be purely local and independent 
of the surrounding; hence it is independent of the spatial dimension of the 
system. Consequently, the on-site interaction U need not be scaled. So we see 
that the scaled Hubbard Hamiltonian 

h = -—= YsY. + U Y1 ™w (i-6) 

has a nontrivial Z — > oo limit, where both terms, the kinetic energy and 
the interaction, are of the same order of magnitude and are thereby able to 
compete; here (i, j) denotes NN sites i and j. It is this competition between the 
two terms which leads to interesting many-body physics. Mathematically this 
is expressed by the fact that the generic matrix elements of the commutator 
between the kinetic and the interaction part of the Hamiltonian do not vanish 
in the d — > oo limit. 

The quantum scaling (1.5) was determined within a fc-space formulation, 
but it can also be derived within a position-space formulation as will be dis- 
cussed next. 



4 In the following we set the Planck constant ft, the Boltzmann constant ks, and 
the lattice constant a equal to unity. 
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1.2.2 Simplifications of the Many-Body Perturbation Theory 

The most important consequence of the scaling (1.5) is the fact that it leads to 
significant simplifications in the investigation of Hubbard-type lattice models 
[12,14-19]. To understand this point better we take a look at the perturbation 
theory in terms of U. At T = and U = the kinetic energy of the electrons 
is given by 

Here = (c^.c J(T )o is the one-particle density matrix which can be inter- 
preted as the probability amplitude for hopping from site j to site i. The 
square of its absolute value is proportional to the probability for an electron 
to hop from j to i, i.e., |<7y CT | 2 ~ 1/Z ~ 1/d, since site j has 0{d) NN sites i. 
The sum of \gij a | 2 over all NN sites i of j must then yield a Z or d independent 
constant. In the limit rf^oowc therefore find 




Since the sum over NN sites in (1.7) is of 0(d), the NN hopping amplitude 
t must obviously be scaled according to (1.5) for E® in to remain finite in the 
limit d,Z — > oo. Hence, as expected, a real-space formulation yields the same 
results for the required scaling of the hopping amplitude. 

The one-particle Green function G^ ;Cr (w) of the non- interacting system 
obeys the same scaling as g\^ a - This follows directly from its definition 

G%^{t) = -{Tc ia {t)c% (0)) , (1.9) 

where T is the time ordering operator, and the time evolution of the operators 
is given by the Hciscnberg representation. The one-particle density matrix is 
obtained as gfj a — lim t ^ - G^- CT (i). If g^j obeys (1.8) the one-particle Green 
function must follow the same scaling at all times since this property does not 
dependent on the time evolution and the quantum mechanical representation. 
The Fourier transform (w) also preserves this property. 

It is important to realize that, although G®j CT ~ l/Vd vanishes for d — > oo, 
the particles are not localized, but are still mobile. Indeed, even in the limit 
d — > oo the off-diagonal elements of C-j contribute, since a particle may hop 

to d nearest neighbors with reduced amplitude t* j\Pld. For general i,j one 
finds [20, 15] 

G° ijta ^o(l/d^- R ^ 2 ), (1.10) 

where ||i?|| = X^n=i ^ s ^ ne length of i? in the so-called "New York metric" 
(also called "taxi cab metric" , since particles only hop along horizontal or 
vertical lines, never along a diagonal). 

It is the property (1.10) which is the origin of all simplifications arising 
in the limit d — > oo. In particular, it implies the collapse of all connected, 
irreducible perturbation theory diagrams in position space [12,14,15]. In gen- 
eral, any two vertices which are connected by more than two separate paths 5 
will collapse onto the same site. In particular, the external vertices of any 



5 Here a "path" is any sequence of lines in a diagram; they are "separate" when 
they have no lines in common. 
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irreducible self-energy diagram are always connected by three separate paths 
and hence always collapse. As a consequence the full irreducible self-energy 
becomes a purely local quantity [12], but retains its dynamics [14] 

Eij,a{u) £ii,<r(u)Sij- (1.11a) 

In the paramagnetic phase we may write Su t(T (cj) = The Fourier trans- 

form of is then momentum-independent 

^(fc, w ) d =°° E a (w). (1.11b) 

This leads to tremendous simplifications in all many-body calculations for the 
Hubbard model and related models. It should be noted that a fc-independence 
of E was sometimes assumed as a convenient approximation ( "local approxi- 
mation") [21-23]. Here we identified the limit where this is indeed exact. 



1.2.3 Interactions Beyond the On-Site Interaction 

In the case of more general interactions than the Hubbard interaction, e.g., 
nearest-neighbor interactions such as 

Hnn = J2J2 Vea'h^, (1.12) 

the interaction constant has to be scaled, too, in the limit d — > oo. In the case 
of (1.12), which has the form of a classical interaction, the "classical" scaling 

y<ja' -+ ^ (1.13) 

is required. Of course, the propagator still has the dependence (1.10). 

Due to (1.13) all contributions, except for the Hartree term, are found 
to vanish in d = oo [14]. Hence nonlocal interactions only contribute via 
their Hartree contribution, which is purely static. This gives the Hubbard 
interaction a unique role: of all interactions for fermionic lattice models only 
the purely local Hubbard interaction remains dynamical in the limit d — > oo 
[14]. 



1.2.4 Single-Particle Propagator 

Due to the fc-independence of the irreducible self-energy, (1.11b), the one- 
particle propagator of an interacting lattice fermion system is given by 

G *." (w) = 7, f a-), r u,Y (L14) 

Most importantly, the k dependence of G k (y>) comes entirely from the energy 
dispersion e k of the non-interacting particles. This means that for a homoge- 
neous system with the propagator 

G ij>a {u) = L- 1 ^ G k ^)e ik - (Ri - Ri) (1.15) 
k 

its local part, Ga tfT , has the form [16] 
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G iiA .) = L -1 £ G kA .) = J dE u _ E N ~W sAu) (1.16a) 

fc — oo 

= G a (w). (1.16b) 

In the following we will limit our discussion to the paramagnetic phase and 
omit the spin index. The spectral function of the interacting system (often 
referred to as the "density of states" (DOS) as in the non-interacting case) is 
given by 

A{to) = - 1 ImG(w + i0+); (1.16c) 

7T 

for U = one has A(uj) = N(co). In the limit d — > oo two quantities then play 
the most important role: the local propagator G(lo) and the self-energy E(uj). 



1.2.5 Consequences of the Momentum Independence of the 
Self-Energy 

We now discuss some more consequences of the fe-independence of the self- 
energy as derived by Muller-Hartmann [16]. Let us consider the Hubbard 
model, or any one of its generalizations, in the paramagnetic phase, i.e., with- 
out a broken symmetry. At T — the one-particle Green function (1.14) then 
reads 

G fc (w) = 1 ^-r. (1.17) 

In general, even when S(oj) is fc-dependent, the Fermi surface is defined by 
the lo — limit of the denominator of (1.17) as 

e k + E k {0) = E F . (1.18a) 

According to Luttinger and Ward [24] the volume within the Fermi surface is 
not changed by interactions, provided the effect of the latter can be treated 
in infinite-order perturbation theory (i.e., no broken symmetry). This is ex- 
pressed by 

n = J2 e \ E F - efc - £fc(0)], (1.18b) 

her 

where n is the particle density and 6{x) is the step function. In general, the 
fc-dependence of S k (0) in (1.18a) implies that, in spite of (1.18b), the shape 
of the Fermi surface of the interacting system will be quite different from 
that of the non-interacting system (except for the fully rotation invariant case 
£fe ~ k 2 ). For lattice fermion models in d < oo, with S k (co) = (1.18a) 
implies that the Fermi surface itself (and hence the volume enclosed) is not 
changed by interactions. 6 The Fermi energy is simply shifted uniformly from 
its non-interacting value E F , i.e., Ep = E F + S(0), to keep n in (1.18b) 
constant. From (1.16a) we thus conclude that the u> = value of the local 
propagator, G(0), and hence of the spectral function, A(0) = — ^lmG(iO + ), 
is not changed by interactions. Rcnormalizations of A(0) can only come from 
a fc-dependence of S, i. e., if dS/dk ^ 0. 

For lo — > the self-energy has the property 



6 In d = oo limit the notion of a Fermi surface of a lattice system is complicated 
by the fact that the dispersion e fe is not a simple smooth function. 
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Im E(oj) oc u 2 (1.18c) 
which implies quasiparticle (Fermi liquid) behavior. The effective mass 



to* _ dS 
m du 



= 1 + 



LJ=0 



i r 

k J-< 



dijj 



lmS(uj + iO 



> 1 



is seen to be enhanced. In particular, the momentum distribution 

rO 



i r 

rife = — / (iw ImGfc(w) 

^* J-oo 



(1.18d) 



(1.19) 



has a discontinuity at the Fermi surface, given by n k - — n k + = (m*/m) 1 , 
± 



where k F = kp ± + . 



1.3 Dynamical Mean-Field Theory for Correlated 
Lattice Fermions 

The limit of high spatial dimensions d or coordination number Z provides 
the basis for the construction of a comprehensive mean-field theory for lattice 
fermions which is diagrammatically controlled and whose free energy has no 
unphysical singularities. It starts from the scaled Hamiltonian (1.6) and makes 
use of the simplifications in the many-body perturbation theory discussed in 
sec. 1.2.2. There we found that the local propagator G(u>), i.e., the probability 
amplitude for an electron to return to a lattice site, and the local, but fully 
dynamical self-energy S{uj) are the most important quantities in such a theory. 
Since the self-energy is a dynamical variable (in contrast to the Hartree-Fock 
theory where it is merely an average, static potential) the resulting mean-field 
theory is also dynamical and can thus describe genuine correlation effects such 
as the Mott-Hubbard metal-insulator transition. 

The self-consistency equations of this dynamical mean-field theory (DMFT) 
for correlated lattice fermions can be derived in different ways. All derivations 
make use of the fact that in the limit of high spatial dimensions Hubbard-type 
models, i.e., lattice models with a local interaction, reduce to a "dynamical 
single-site problem", where the d-dimensional lattice model is effectively de- 
scribed by the dynamics of the correlated fermions on a single site embed- 
ded in a "bath" provided by the other particles. In particular, the derivation 
by Janis [25,26] is a generalization of the coherent potential approximation 
(CPA) for disordered systems 7 to the Hubbard model. In the following we 
will present today's standard derivation by Georges and Kotliar [31] which 
is based on the mapping of the lattice problem onto a self-consistent single- 
impurity Anderson model; this approach was also employed by Jarrell [32]. 
Although the DMFT equations derived within the CPA approach and the self- 
consistent single-impurity approach, respectively, are identical it is the latter 
formulation which was immediately adopted by the community since it makes 
contact with the theory of quantum impurities and Kondo problems; for a 
review see [33]. This is a well- understood branch of many-body physics [34] 
for whose solution efficient numerical codes had been developed already in the 
1980's, in particular by use of the quantum Monte-Carlo (QMC) method [35]. 



7 The CPA is the best single-site approximation for disordered, non-interacting 
lattice electrons [27-29]; it becomes exact in the limit d,Z^oo [30,8]. 
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d=3, Z=12 




Fig. 1.1. Already in a d — 3 the coordination number Z can be quite large, as in the 
case of a face-centered cubic lattice with Z — 12. In the limit d — > oo, i.e., Z — > oo, 
the many-body lattice problem reduces to that of a single lattice site embedded in 
a dynamical mean field. As shown in the inset electrons can hop onto and off that 
site and interact as in the finite-dimensional Hubbard model. Therefore the DMFT 
describes the dynamics of the interacting electrons correctly. 



1.3.1 Construction of the DMFT as a Self-Consistent 
Single-Impurity Anderson Model 

Following the presentation of Georges, Kotliar, Krauth, and Rozenberg [33] 
the DMFT equations will now be derived using the so-called cavity method 8 . 
This derivation starts by removing one lattice site together with its bonds from 
the rest of the lattice. The remaining lattice, which now contains a cavity, is 
replaced by a particle bath which plays the role of the dynamical mean field 
(see Fig. 1.1). Now comes a physically motivated step: the bath is coupled, 
via a hybridization, to the cavity. The resulting problem then amounts to the 
solution of an effective single-impurity Anderson model where the degrees of 
freedom of the bath, represented by an appropriate hybridization function, 
have to be determined self-consistently. 

To be specific, we consider the partition function in the grand canonical 
ensemble 

Z = f\{Dc* a Dc. laG M-S]. (1-20) 
The action S for the Hubbard model is given by 



S= I dr 







X, c **( r )(a= -M)cio-(r) + 2j%c* CT (T)cjv(r) 



lj a 



(1.21) 



where we use Grassmann variables c* al Ci a . We split the action S into three 
parts 

S = S + AS + S {0 \ (1.22) 
where So is the part containing only variables on site 



We note that the sign of the hopping amplitude tij used here (see the definition 
in (1.1b)) is opposite to that in ref. [33]. 
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S = I dr 

10 



H C CT ( T )(g7 _ M)C0 CT (T) + Uc* ot {T)c 01 -(T)c* Ql (T)col(T) 



(1.23) 

AS contains the hoppings between site and other sites of the lattice i ^ 
AS = / dTy2[tioct(r)c ^T)+t 0i c* 0<T (T)c ia (T)], (1.24) 

and the rest, denoted by S^°\ is the part of the action where the site and 
its bonds are removed, i.e., for i, j ^ one has 



Jo 



dr 



d 



+ U^C* t (T)c it (T)c* i (T)c il (T) 



(1.25) 



We now rewrite the partition function Z as 
Z = I Y\DcQ a Dc 0a cxp[-So} 

x f J| Dc* irT Dc ia exp[-5 (0) ] exp[-AS] (1.26) 

and use the ensemble average 

W(0) = i / II DclD Ci „XeM-SM] (1.27) 

taken with respect to S^ (the action where the site i — is excluded), with 
z(°) being the corresponding partition function. Then the partition function 
roads 

Z = Z<°) fl[Dc* „Dco a exp[-S }(exp[-AS}} {0) . (1.28) 

In the next step we expand the second exponent with respect to the action 
AS. A non-trivial limit d — > oo is obtained by scaling the hopping amplitudes 
tij as described in sec. 1.2.2. Consequently, in the Z — >• oo limit only the 
contribution G^ CT , where 



Gf^i - r 2 ) = - {T T Cj a (t i)c* ka (r 2 )) (o) , 



(1.29) 



or disconnected contributions made of products of G^.'s remain. Applying 
the linked-cluster theorem and collecting only connected contributions in the 
exponential function one obtains the local action 
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S\oc — 



rP Q 

rP 
Jo 

fP rP 

+ dn dr 2 V V tj t% Gf^(n - r 2 )c* 0(T (T 1 )c 0a (T2) 



(1.30) 



where the rescaled hoppings are denoted with a star. Introducing the hy- 
bridization function 



A*(n - 75) = - £ fnf jQ G§l(n r 2 ), 



(1.31) 



and employing the free ("Weiss") mean-field propagator Q a one can express 
the DMFT local action in the following form (here the site index i = is 
omitted for simplicity) 

r P r P 

S loc = - / dn dr 2 V cKr^g- 1 ^ - T 2 )c a (T 2 ) 
Jo Jo a 

rP 

+ U / dT4(T)ct(r)c|(r)c 4 (r), (1.32) 
Jo 



where 



G a 1 {n ~t 2 ) = - (J^- - f?j S TlT2 - A a {n - t 2 ). 



(1.33) 



Finally we need the relation between the Green function G^1(t — r') 
where the site i = is removed and the full lattice Green function, i.e., 



ijcr ~ u ij<? ^iOff^oOcr^Ojcrj 



(1.34) 



which holds for a general lattice. 

In order to obtain the full solution of the lattice problem it is convenient 
to express the relation between the local Green function Gooa = G a and the 
"Weiss" mean field 9 Q' 1 in the form of a Dyson equation 



[G^iWn)}- 1 = [GcriiUn)}- 1 - E a (iu n ) 

= iw n +(J,- A a (iuj n ) - £ a (iLj n ). 
Then the lattice Green function (in fe-space) Gkcr(iu n ) is given by 

G ka {iuj n ) 



iuj n - efe + H - S a (iiJ n ) 



(1.35) 
(1.36) 



(1.37) 



After performing the so-called lattice Hilbert transform we recover the local 
Green function 



In principle, any one of the local functions Ga(ioj n ), S a {iu> n ), or Z\ CT (icj n ) can be 
viewed as a "dynamical mean field" acting on particles on a site, since they all 
appear in the bilinear term of the local action (1.32). 
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G a (iu n ) = VG feCT (iO= J2~. — ^ ^y-— (1.38) 

OO 

= [de (1 . 39 ) 

— oo 

The ionic lattice on which the electrons move, and its structure, are seen 
to enter only via the DOS of the non-interacting electrons. After analytic 
continuation to real frequencies the local ( "fc averaged" ) propagator reads 

Gfe.M = — ^ (1.40) 

This completes the derivation of the self-consistent DMFT equations. 
Namely, the functional integral determining the local propagator 

G a (iuj n ) = -— fjJ_DclDc a [c a (iu n )cl(iu} n )]ejq?[-Si oc ], (1.41) 

where the partition function Z and the local action Si oc are given by 

Z= fl[Dc* a Dc a e X p[-Sio C ] (1-42) 

and (1.32), respectively, together with the expression (1.39) for the lattice 
Green function provide a closed set of equations for the local propagator 
G a (iuj n ) and the self-energy E a (iuj n ). These equations can be solved itera- 
tively: Starting with an initial value for the self-energy one obtains the local 
propagator from (1.39) and thereby the bath Green function Q a (iiv n ) from 
(1.35). This determines the local action (1.32) which is used to solve the 
single-impurity problem (1.41), leading to a new value for the local propaga- 
tor and, by employing the old self-energy, a new bath Green function, etc. The 
single-impurity problem is still a complicated many-body interacting problem 
which cannot, in general, be solved exactly. 



1.3.2 Solution of the Self- Consistency Equations of the DMFT 

Due to the purely on-site interaction in the local action (1.32) the dynamics 
of the Hubbard model, (1.1), remains complicated even in the limit d — > oo. 
Exact evaluations are only feasible when there is no coupling between the 
frequencies as, for example, in the Falicov-Kimball model [36]. This model was 
solved analytically by Brandt and Mielsch [19] soon after the introduction of 
the d -> oo limit [12]. 

In general, the local action (1.32) is the most complicated part of the 
DMFT equations. To solve the self-consistency equations different techniques 
( "impurity solvers" ) have been developed which are either fully numerical and 
"numerically exact" , or semi-analytic and approximate. The numerical solvers 
can be divided into renormalization group techniques such as the numerical 
renormalization group (NRG) [37,38] and the density-matrix renormalization 
group (DMRG) [39], exact diagonalization (ED) [40-42], and methods based 
on the stochastic sampling of quantum and thermal averages, i.e., quantum 
Monte-Carlo (QMC) techniques such as the Hirsch-Fye QMC algorithm [32, 
43,44,33] and continuous-time (CT) QMC [45-47]. 
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Semi-analytic approximations such as the iterated perturbation theory 
(IPT) [31,48,33], the non-crossing approximation (NCA) [33,49], the fluc- 
tuation exchange approximation (FLEX) [50-53] , the local moment approach 
(LMA) [54, 55] , and the parquet approximation [56] can also provide valuable 
insight. 

It quickly turned out that the DMFT is a powerful tool for the investigation 
of electronic systems with strong correlations. It provides a non-perturbative 
and thermodynamically consistent approximation scheme for fmitc-dimension- 
al systems which is particularly valuable for the study of intermediate-coupling 
problems where perturbative techniques fail [57,33,58,59,8]. 

1.4 The Mott-Hubbard Metal-Insulator Transition 

The correlation induced transition between a paramagnetic metal and a 
paramagnetic insulator, referred to as "Mott-Hubbard metal-insulator tran- 
sition (MIT)", is one of the most intriguing phenomena in condensed matter 
physics [60-62]. This transition is a consequence of the competition between 
the kinetic energy of the electrons and their local interaction U. Namely, the 
kinetic energy prefers the electrons to move (a wave effect) which leads to 
doubly occupied sites and thereby to interactions between the electrons (a 
particle effect). For large values of U the doubly occupied sites become ener- 
getically very costly. The system may reduce its total energy by localizing the 
electrons. Hence the Mott transition is a localization-delocalization transition, 
demonstrating the particle- wave duality of electrons [59] . 

Mott-Hubbard MITs are found, for example, in transition metal oxides 
with partially filled bands. For such systems band theory typically predicts 
metallic behavior. The most famous example is V2O3 doped with Cr [63-65]. 
In particular, below T = 380 K the metal-insulator transition in paramagnetic 
(Vo.geCro. 04)203 is of first order [64], with discontinuities in the lattice param- 
eters and in the conductivity. However, the two phases remain isostructural. 

Making use of the half-filled, single-band Hubbard model (1.1) the Mott- 
Hubbard MIT was studied intensively in the past [5,65,60-62]. Important 
early results were obtained by Hubbard [66] within a Green function decou- 
pling scheme, and by Brinkman and Rice [67] within the Gutzwiller variational 
method [4], both at T = 0. Hubbard's approach yields a continuous splitting 
of the band into a lower and upper Hubbard band, but cannot describe quasi- 
particle features. By contrast, the Gutzwiller-Brinkman-Rice approach gives 
a good description of the low-energy, quasiparticle behavior, but cannot re- 
produce the upper and lower Hubbard bands. In the latter approach the MIT 
is signalled by the disappearance of the quasiparticle peak. 

To solve this problem the DMFT has been extremely valuable since it 
provided detailed insights into the nature of the Mott-Hubbard MIT for all 
values of the interaction U and temperature T [33,59]. 

1.4.1 DMFT and the Three-Peak Structure of the Spectral 
Function 

The Mott-Hubbard MIT is monitored by the spectral function A(oj) = 
— -ImC?(u> + i0 + ) of the correlated electrons 10 [33]. While at small U the 



In the following we only consider the paramagnetic phase, whereas magnetic order 
is assumed to be suppressed ("frustrated"). 
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system can be described by coherent quasiparticles whose DOS still resembles 
that of the free electrons, the spectrum in the Mott insulator state consists 
of two separate incoherent "Hubbard bands" whose centers are separated 
approximately by the energy U. The latter originate from atomic- like excita- 
tions at the energies ±C//2 broadened by the hopping of electrons away from 
the atom. At intermediate values of U the spectrum then has a characteristic 
three-peak structure as in the single-impurity Anderson model, which includes 
both the atomic features (i.e., Hubbard bands) and the narrow quasiparticle 
peak at low excitation energies, near uj = 0. This corresponds to a strongly 
correlated metal. The general structure of the spectrum (lower Hubbard band, 
quasiparticle peak, upper Hubbard band) is rather insensitive to the specific 
form of the DOS of the non-interacting electrons. The width of the quasipar- 
ticle peak vanishes for U — > U c2 (T). On decreasing U, the transition from the 
insulator to the metal occurs at a lower critical value U c i(T), where the gap 
vanishes. 

It is important to note that the three-peak spectrum originates from a 
lattice model with only one type of electrons. This is in contrast to the single- 
impurity Anderson model whose spectrum shows very similar features, but is 
due to two types of electrons — the localized orbital at the impurity site and 
the free conduction band. Therefore the screening of the magnetic moment 
which gives rise to the Kondo effect in impurity systems has a somewhat 
different origin in interacting lattice systems. Namely, as explained by the 
DMFT the electrons provide both the local moments and the electrons which 
screen these moments [33]. 

Interestingly, for any typical spectral function A(u) with three peaks, 
Kramers-Kronig relations and the DMFT self-consistency equations imply 
that the self-energy E{u>) abruptly changes slope inside the central peak at 
some frequency cj* [68] , once at positive and once at negative frequency. While 
this behavior is not visible in A(u>) itself, it leads to "kinks" in the effec- 
tive dispersion relation Ek of one-particle excitations, which is defined as 
the frequency for which the momentum-resolved spectral function A(k,uj) = 
— ImG(fc, uj)/tt = — (l/7r) lm[l/((jj+fi—e-k — S(oj))] is maximal. For frequencies 
below lo+ the dispersion is given by Fermi-liquid (FL) theory, Ek = ^FLCfe, 
where Z FL = (dBjeS(u)/duj) u= Q is the FL renormalization parameter. The 
FL regime terminates at the kink energy scale w*. This energy cannot be 
obtained within FL theory itself. Namely, it is determined by Zfl and the 
non-interacting DOS, e.g., = OAlZpi^D, where D is an energy scale of 
the non- interacting system such as half the bandwidth [68]. Above w+ the 
dispersion is given by a different renormalization with a small offset, Ek = 
Zcptk + const, where Zcp is the weight of the central peak of A(co). This the- 
ory explains kinks in the slope of the dispersion as a direct consequence of the 
electronic interaction [68] . The same mechanism may also lead to kinks in the 
low-temperature electronic specific heat [69] . These kinks have also been linked 
to maxima in the spin susceptibility [70]. Of course, additional kinks in the 
electronic dispersion may also arise from the coupling of electrons to bosonic 
degrees of freedom, such as phonons [71,72] or spin fluctuations [73,74]. In- 
terestingly, recent experiments [75] have found evidence for kinks in Ni(110), 
which may be due to the electronic mechanism discussed here. 

The evolution of the spectral function of the half-filled frustrated Hubbard 
model at a finite temperature is shown in Fig. 1.2. This temperature is above 
the temperature of the critical point so that there is no real transition but only 
a crossover from a metallic-like to an insulating-like solution. The height of the 
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Fig. 1.2. Spectral function for the half-filled Hubbard model for various values of 
U at T = 0.0276 W (W: bandwidth) in the crossover region. The crossover from 
the metal to the insulator occurs via a gradual suppression of the quasiparticle peak 
at cj = 0. The inset shows the U dependence of A(uj = 0), in particular the rapid 
decrease for U « 1.1 W; from Ref. [76]. 



quasiparticle peak at the Fermi energy is no longer fixed at its zero temper- 
ature value. This is due to the temperature dependent imaginary part of the 
self-energy. The spectral weight of the quasiparticle peak is seen to be gradu- 
ally redistributed and shifted to the upper (lower) edge of the lower (upper) 
Hubbard band. The inset of Fig. 1.2 shows the [/-dependence of the spectral 
function at zero frequency A(w = Q). For higher values of U the spectral den- 
sity at the Fermi level is still finite and vanishes only in the limit U — > oo. 
For the insulating phase DMFT predicts the filling of the Mott-Hubbard gap 
with increasing temperature. This is due to the fact that the insulator and 
the metal are not distinct phases in the crossover regime, implying that the 
insulator has a finite spectral weight at the Fermi level. This behavior has 
been detected experimentally by photoemission experiments [77]. 

Altogether, the thermodynamic transition line U C (T) corresponding to the 
Mott-Hubbard MIT is found to be of first order at finite temperatures, and 
is associated with a hysteresis region in the interaction range U c \ < U < 
U C 2, where U c i and C/ c2 are the interaction values at which the insulating 
and metallic solution, respectively, vanish [33,37,78,79,76,80]. As shown in 
Fig. 1.3 the hysteresis region terminates at a critical end point. At higher 
temperatures the transition changes into a smooth crossover from a bad metal 
to a bad insulator. At half filling and for bipartite lattices in dimensions 
d > 2 (in d = 2 only at T = 0) the paramagnetic phase is, in fact, unstable 
against antiferromagnetic long-range order. The metal-insulator transition is 
then completely hidden by the antiferromagnetic insulating phase [81]. 

In Fig. 1.3 it is seen that the slope of the phase transition line U c is negative 
down to T = 0, which implies that for constant interaction U the metallic 
phase can be reached from the insulator by decreasing the temperature T, i.e., 
by cooling. This anomalous behavior (which corresponds to the Pomeranchuk 
effect [82] in 3 He, if we associate solid 3 He with the insulator and liquid 3 Hc 
with the metal) can be understood from the Clausius-Clapeyron equation 
dll/dT — AS/ AD. Here AS is the difference between the entropy in the metal 
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Fig. 1.3. Phase diagram of the Mott-Hubbard MIT showing the metallic phase and 
the insulating phase, respectively, at temperatures below the critical end point, as 
well as a coexistence region; from Ref. [80]. 

and in the insulator, and AD is the difference between the number of doubly 
occupied sites in the two phases. Within the DMFT there is no exchange 
coupling J between the spins of the electrons in the insulator, since the scaling 
(1.5) implies J oc —t 2 /U oc 1/d — > for d — > oo. Hence the insulating state is 
macroscopically degenerate, with entropy Si ns = &b In 2 per electron down to 
T = 0. This is larger than the entropy S mc t oc T per electron in the Landau 
Fermi liquid describing the metal, i.e., AS = S mc t — S- ms < 0. At the same 
time the number of doubly occupied sites is lower in the insulator than in the 
metal, i.e., AD = D met — D- lns > 0. The Clausius-Clapeyron equation then 
implies that the phase-transition line T vs. U has indeed a negative slope down 
to T = 0. However, this is an artifact of the DMFT. Namely, there will usually 
exist an exchange coupling between the electrons which leads to a vanishing 
entropy of the insulator for T — > 0. Since the entropy of the insulator vanishes 
faster than linearly with the temperature, the difference AS = S mct — Si ns 
eventually becomes positive, whereby the slope also becomes positive 11 ; this 
is indeed observed in cluster DMFT calculations [83] . However, since AS — > 
for T — > the phase boundary must eventually terminate at T = with 
infinite slope. 

1.5 Theory of Electronic Correlations in Materials 
1.5.1 The LDA+DMFT Approach 

Although the Hubbard model is able to explain basic features of the phase 
diagram of correlated electrons it cannot describe the physics of real materials 
in any detail. Clearly, realistic theories must take into account the explicit 
electronic and lattice structure of the systems. 

Until recently the electronic properties of solids were investigated by two 
essentially separate communities, one using model Hamiltonians in conjunc- 
tion with many-body techniques, the other employing density functional the- 
ory (DFT) [84,85]. DFT and its local density approximation (LDA) have 

11 Here we assume for simplicity that the metal remains a Fermi liquid and the 
insulator stays paramagnetic down to the lowest temperatures. 
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the advantage of being ab initio approaches which do not require empirical 
parameters as input. Indeed, they are highly successful techniques for the cal- 
culation of the electronic structure of real materials [86]. However, in practice 
DFT/LDA is seriously restricted in its ability to describe strongly correlated 
materials where the on-site Coulomb interaction is comparable with the band 
width. Here, the model Hamiltonian approach is more powerful since there 
exist systematic theoretical techniques to investigate the many-electron prob- 
lem with increasing accuracy. Nevertheless, the uncertainty in the choice of 
the model parameters and the technical complexity of the correlation prob- 
lem itself prevent the model Hamiltonian approach from being a suitable tool 
for studying real materials. The two approaches are therefore complementary. 
In view of the individual power of DFT/LDA and the model Hamiltonian 
approach, respectively, a combination of these techniques for ab initio investi- 
gations of correlated materials including, for example, /-electron systems and 
Mott insulators, would be highly desirable. One of the first successful attempts 
in this direction was the LDA+ U method [87, 88] , which combines LDA with 
a static, i.e., Hartree-Fock-like, mean- field approximation for a multi-band 
Anderson lattice model with interacting and non-interacting orbitals. This 
method proved to be a very useful tool in the study of long-range ordered, 
insulating states of transition metals and rare-earth compounds. However, 
the paramagnetic metallic phase of correlated electron systems such as high- 
temperature superconductors and heavy-fermion systems clearly requires a 
treatment that goes beyond a static mean-field approximation and includes 
dynamical effects, i.e., the frequency dependence of the self-energy. 

Here the recently developed LDA+DMFT method, a new computational 
scheme which merges electronic band structure calculations and the DMFT, 
has proved to be a breakthrough [89-94,49,59,52,95-98]. Starting from con- 
ventional band structure calculations in the LDA the correlations are taken 
into account by the Hubbard interaction and a Hund's rule coupling term. 
The resulting DMFT equations are solved numerically, e.g., with a quantum 
Monte-Carlo (QMC) algorithm. By construction, LDA+DMFT includes the 
correct quasiparticlc physics and the corresponding energetics. It also repro- 
duces the LDA results in the limit of weak Coulomb interaction U. More im- 
portantly, LDA+DMFT correctly describes the correlation induced dynamics 
near a Mott-Hubbard MIT and beyond. Thus, LDA+DMFT is able to account 
for the physics at all values of the Coulomb interaction and doping level. 

In the LDA+DMFT approach the LDA band structure is expressed by 
a one-particle Hamiltonian H® BA , and is then supplemented by the local 
Coulomb repulsion U and Hund's rule exchange J. This leads to a material 
specific generalization of the one-band model Hamiltonian 

m i i,m^m' , cr.tr ' 

, . (L43) 

Here m and ml enumerate those orbitals for which the interaction be- 
tween the electrons is explicitly included, e.g., the three ti g orbitals of the 3d 
electrons of transition metal ions or the 4/ orbitals in the case of rare earth 
elements. The interaction parameters are related by V = U — 2 J which holds 
exactly for degenerate orbitals and is a good approximation for t 2g electrons. 
The actual values for U and V can be calculated by constrained LDA [49] . 

In the one-particle part of the Hamiltonian 
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^LDA = #LDA - Yl S ACd (1 ' 44) 

i ma 

the energy term containing Aej is a shift of the one-particle potential of the 
interacting orbitals. It cancels the Coulomb contribution to the LDA results, 
and can also be calculated by constrained LDA [49]. 

Within the LDA+DMFT scheme the self-consistency condition connecting 
the self-energy E and the Green function G at frequency u> reads: 

G qm>q ' m >(w) =77- / ^([wl+Ml-^VW-^M]" 1 ) (1-45) 

Here, 1 is the unit matrix, /i the chemical potential, i?L DA (fc) is the orbital 
matrix of the LDA Hamiltonian derived, for example, in a linearized muffin-tin 
orbital (LMTO) basis, S(w) denotes the self-energy matrix which is nonzero 
only between the interacting orbitals, and [...] _1 implies the inversion of the 
matrix with elements n (=qm), n! (—q'm 1 ), where q and m are the indices of the 
atom in the primitive cell and of the orbital, respectively 12 . The integration 
extends over the Brillouin zone with volume Vb- 

For cubic transition metal oxides Eq. (1.45) can be simplified to 

GM= ^-27 M) = | & — 0- (1.46) 

provided the degenerate t 2g orbitals crossing the Fermi level are well separated 
from the other orbitals [49]. For non-cubic systems the degeneracy is lifted. 
In this case eq. (1.46) is an approximation where different E m (u), N^(e) and 
G m (uj) have to be used for the three non-degenerate t 2g orbitals. 

The Hamiltonian (1.43) is diagonalizcd within the DMFT where, for ex- 
ample, quantum Monte-Carlo (QMC) techniques [35] are used to solve the 
self-consistency equations. From the imaginary time QMC Green function we 
calculate the physical (real frequency) spectral function with the maximum 
entropy method [99]. 

During the last few years the LDA+DMFT and other DMFT based com- 
putational schemes have provided great progress in the understanding of the 
electronic, magnetic and structural properties of many correlated electron ma- 
terials. These materials range from 3d transition metals and their oxides, and 
/ electron systems, all the way to Heusler alloys, ferromagnetic half-metals, 
fullerenes, and zeolites [52,95-98]. Nevertheless, this framework still needs to 
be considerably improved before it becomes a truly comprehensive ab initio 
approach for complex correlated matter with predictive power. In particular, 
the interface between the band-structure and the many-body components of 
the approach needs to be optimized. This includes, for example, a solution 
of the double counting correction problem and a fully self-consistent treat- 
ment of the spin, orbital, and charge densities. Another important goal are 
realistic computations of free energies and forces, and the development of 
efficient methods to treat non-local correlations with quantum cluster meth- 
ods [100-102]. 

1.5.2 Single-Particle Spectrum of Correlated Electrons in 
Materials 

Transition metal oxides are an ideal laboratory for the study of electronic cor- 
relations in solids. Among these materials, cubic perovskites have the simplest 

12 We note that H^oa may include additional non-interacting orbitals. 
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crystal structure and thus may be viewed as a starting point for understand- 
ing the electronic properties of more complex systems. Typically, the 3d states 
in those materials form comparatively narrow bands with width W~ 2—3 eV, 
which leads to strong Coulomb correlations between the electrons. Particu- 
larly simple are transition metal oxides with a 3d 1 configuration since, among 
others, they do not show a complicated multiplet structure. 

Photocmission spectra provide a direct experimental tool to study the elec- 
tronic structure and spectral properties of electronically correlated materials. 
Intensive experimental investigations of spectral and transport properties of 
strongly correlated 3d 1 transition metal oxides started with investigations by 
Fujimori et al. [103]. These authors observed a pronounced lower Hubbard 
band in the photoemission spectra which cannot be explained by conven- 
tional band structure theory. In photoemission spectroscopy (PES) a photon 
of a given energy is used to emit an electron whose properties (energy, angu- 
lar distribution) are measured by a detector. Angular resolved PES is referred 
to as ARPES. These techniques measure the occupied electronic states, i.e., 
those states which are described by the full spectral function multiplied by 
the Fermi function f(uj,T). By contrast, inverse photoemission spectroscopy 
(IPES) measures the unoccupied electronic states, i.e., the states described 
by the full spectral function of a material multiplied by 1 — f(ui,T). IPES is 
harder to perform and not as accurate as PES. But in many situations in- 
formation about the unoccupied states is also available by X-ray absorption 
spectroscopy (XAS). 

Spectroscopic techniques provide very valuable information about corre- 
lated electronic systems since they can directly measure the spectral function 
of a material, a quantity which can also be directly calculated as discussed 
in sec. 1.4. In particular, photoemission techniques allow one to detect the 
correlation induced shift of spectral weight . In the following we will illustrate 
the computation of the fc-integrated electronic spectra of correlated materials 
within the LDA+DMFT scheme by investigating the two simple transition 
metal oxides SrVC-3 and CaVC>3. 

Sr,C ai _ x V0 3 

SrV03 and CaV03 are simple transition metal compounds with a 3d 1 config- 
uration. The main effect of the substitution of Sr ions by the isovalent, but 
smaller, Ca ions is to decrease the V-O-V angle from 8 = 180° in SrVC-3 to 
8 w 162° in the orthorhombically distorted structure of CaV03. Remarkably 
this rather strong bond bending results only in a 4% decrease of the one- 
particle bandwidth W and thus in a correspondingly small increase of the 
ratio U/W as one moves from SrV0 3 to CaV0 3 [104,105]. 

LDA+DMFT(QMC) spectral functions of SrV0 3 and CaV0 3 were calcu- 
lated by Sekiyama et al. [104] by starting from the respective LDA DOS of 
the two materials; they are shown in Fig. 1.4. These spectra show genuine 
correlation effects, i.e., the formation of lower Hubbard bands at about 1.5 eV 
and upper Hubbard bands at about 2.5 eV, with well-pronounced quasiparti- 
cle peaks at the Fermi energy. Therefore both SrVC>3 and CaV03 are strongly 
correlated metals. The small difference of the LDA bandwidth of SrV03 and 
CaV03 is only reflected in some additional transfer of spectral weight from 
the quasiparticle peak to the Hubbard bands, and minor differences in the po- 
sitions of the Hubbard bands. The DOS of the two systems shown in Fig. 1.4 
arc quite similar. In fact, SrV0 3 is slightly less correlated than CaV0 3 , in 
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Fig. 1.4. LDA+DMFT(QMC) spectral functions of SrV0 3 (solid line) and CaVOa 
(dashed line) calculated at T=300 K; inset: effect of temperature in the case of 
CaVOs; after Ref. [104]. 



accord with their different LDA bandwidths. The inset of Fig. 1.4 shows that 
the effect of temperature on the spectrum is weak for T < 700 K. Detailed 
spectra of SrVC>3 and CaVC>3 were also computed by Pavarini et al. [106]. 

Since the three tig orbitals of this simple 3d 1 material are almost degen- 
erate the spectral function has the same three-peak structure as that of the 
one-band Hubbard model shown in Fig. 1.2. The temperature induced de- 
crease of the quasiparticle peak height is also clearly seen. We note that the 
actual form of the spectrum no longer resembles the LDA DOS used as input, 
i.e., it essentially depends only on the first three energy moments of the LDA 
DOS (electron density, average energy, band width). 

In the left panel of Fig. 1.5 the LDA+DMFT(QMC) spectral functions 
at 300K are compared with experimental high-resolution bulk PES. For this 
purpose the full theoretical spectra were multiplied with the Fermi function 
at the experimental temperature (20 K) and were Gauss broadened with the 
experimental resolution of 0.1 eV [104]. The quasiparticle peaks in theory and 
experiment are seen to be in very good agreement. In particular, their height 
and width are almost identical for both SrV03 and CaV03. The difference in 
the positions of the lower Hubbard bands may be partly due to (i) the sub- 
traction of the estimated oxygen contribution (which may also remove some 
3d spectral weight below —2 eV), and (ii) uncertainties in the ah initio calcu- 
lation of the local Coulomb interaction strength. In the right panel of Fig. 1.5 
comparison is made with XAS data of Inoue et al. [107]. To this end the 
full LDA+DMFT spectrum was multiplied with the inverse Fermi function 
at 80K and was then Gauss broadened with the experimental resolution of 
0.36 eV [108]. The overall agreement of the weights and positions of the quasi- 
particle and upper ti g Hubbard band is good, including the tendencies when 
going from SrV03 to CaVOs (in fact, Cao.9Sro.1VO3 in the experiment). For 
CaV03 the weight of the quasiparticle peak is somewhat lower than in the 
experiment. In contrast to one-band Hubbard model calculations, the mate- 
rial specific results reproduce the strong asymmetry around the Fermi energy 
w.r.t. weights and bandwidths. 
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Fig. 1.5. Comparison of the calculated, parameter-free LDA+DMFT(QMC) spec- 
tral functions for SrVOa (solid line) and CaV03 (dashed line) with experiment. Left: 
Bulk-sensitive high-resolution PES (SrVOg: circles; CaVC>3: rectangles). Right: Is 
XAS for SrVC>3 (diamonds) and Cao.9Sro.1VO3 (triangles) [107]. Horizontal line: 
experimental subtraction of the background intensity; after Ref. [105]. 



The experimentally determined spectra of SrVC-3 and CaV03 and the 
good agreement with parameter-free LDA+DMFT calculations confirm the 
existence of a pronounced three-peak structure in a correlated bulk material. 
Although the DMFT had predicted such a behavior for the Hubbard model 
[33] it was not clear whether the DMFT would really be able to describe real 
materials. Now it is clear that the three-peak structure not only occurs in 
single-impurity Anderson models or the DMFT for the Hubbard model, but 
is also a characteristic feature of correlated bulk matter in d = 3. 



1.6 Electronic Correlations and Disorder 

The properties of real materials are influenced not only by the interaction be- 
tween the electrons in the periodic crystal lattice, but also by the presence of 
randomness, e.g., impurities and lattice defects [109]. In particular, Coulomb 
correlations and disorder are separately driving forces behind metal-insulator 
transitions (MITs) connected with the localization and derealization of par- 
ticles. The Mott-Hubbard MIT is caused by the electronic repulsion [110,5,61] 
and is characterized by the opening of a gap in the density of states at the 
Fermi level. By contrast, the Anderson localization transition is due to co- 
herent backscattering of non-interacting particles from randomly distributed 
impurities [111, 109]. At the Anderson transition the character of the spectrum 
at the Fermi level changes from a continuous to a dense point spectrum. Both 
MITs can be characterized by a single quantity, the local density of states 
(LDOS). Although the LDOS is not an order parameter associated with a 
symmetry-breaking phase transition, it can distinguish between a metal and 
an insulator. 
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The DMFT can easily be extended to study correlated lattice electrons 
with local disorder [112-117,8]. For this purpose a singlc-particlc term with 
random local energies ej is added to the Hubbard model, leading to the 
Anderson-Hubbard model 



The ionic energies ej describe the local, quenched disorder acting on the motion 
of the electrons. They are drawn from a probability distribution function 
Vfa), which can be a continuous or a multi-modal function. 

The DMFT provides a valuable, non-perturbative theoretical framework 
also for the investigation of correlated electrons in the presence of disorder. 
If in the DMFT the effect of local disorder is taken into account through 
the arithmetic mean of the LDOS one obtains, in the absence of interactions 
(U = 0), the coherent potential approximation (CPA) [30,118], which does 
not describe the physics of Anderson localization. To overcome this deficiency 
Dobrosavljevic and collaborators formulated a variant of the DMFT where 
the geometrically averaged LDOS is computed from the solutions of the self- 
consistent stochastic DMFT equations [119] which is then incorporated into 
the self-consistency cycle [120]. Thereby a mean-field theory of Anderson local- 
ization can be derived which reproduces many of the expected features of the 
disorder-driven MIT for non-interacting electrons [120]. This scheme uses only 
one-particle quantities and can therefore easily be included in the DMFT to 
treat disordered electrons in the presence of phonons [121] or Coulomb corre- 
lations [122,123]. In particular, the DMFT with geometrical averaging allows 
one to compute the phase diagram for the Anderson-Hubbard model [122,123]. 
For a continuous disorder distribution function and a half filled band one finds 
that the metallic phase is enhanced at small and intermediate values of the 
interaction and disorder, but that the metallicity is eventually destroyed when 
the disorder is strong enough [122]. Surprisingly, the Mott and Anderson in- 
sulators are continuously connected. 

New interesting phenomena are also expected in correlated electron sys- 
tems when the disorder distribution function has a binary-alloy form. Namely, 
it was predicted that a disorder induced splitting of the band ("alloy-band 
splitting" ) can enhance the critical temperature for the onset of itinerant fer- 
romagnctism [124,125]. Another direct consequence of the alloy-band splitting 
is the fact that, if the alloy subband is effectively half-filled, a Mott-Hubbard 
metal-insulator transition can occur at non-integer filling [126]. The spectral- 
weight transfer in correlated systems with binary-alloy disorder was also in- 
vestigated away from the metal-insulator transition regime [127,128]. For the 
periodic Anderson model with the binary-alloy disorder analogous behavior 
was predicted [129]. 

1.7 DMFT for Correlated Bosons in Optical Lattices 

The observation of Bose-Einstein condensation (BEC) in ultra-cold atomic 
gases has greatly stimulated research into the properties of this fascinat- 
ing quantum state of matter [130]. In particular, experiments with alkali 
atoms confined in optical lattices [131-133] have renewed the theoretical in- 
terest [134-136] in the physics of strongly correlated bosons on lattices, which 




(1.47) 
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promises significant new insights and even applications in fields such as quan- 
tum computing [137]. 

A lattice model of interacting bosons with purely local interaction, the 
bosonic Hubbard model, has the form 

# BH - J2 ^Mh + \ E - !)' ( L48 ) 

(id) R. 

where hi and b\ are bosonic operators. In the ground state two characteris- 
tically different phases are expected to occur: a bosonic incompressible Mott 
phase with a correlation gap and a compressible superfluid phase charac- 
terized by a non-vanishing expectation value which serves as an order 
parameter [138]. 

The construction of a DMFT for the bosonic Hubbard model which — as in 
the case of the fermionic DMFT — becomes exact in the limit d or Z — > oo and 
is valid at all temperatures is made complicated by the fact that the system 
can Bose condense below a Bose- Einstein condensation temperature Tbec- 
This has immediate consequences for the expectation value of the kinetic 
energy in (1.48). Namely, while in the normal state only the product b\bj has 
a finite expectation value, in the Bose condensed phase even (bj) (the order 
parameter) is non-zero. For the expectation value of the kinetic energy in 
(1.48) to remain finite in the limit Z — > oo the hopping amplitude then has 
to be scaled differently in the two phases. In the normal phase the hopping 
amplitude needs to be scaled as in the fermionic case ("quantum scaling"), i.e., 
as tij — t*j/\/Z, while in the condensed phase a classical scaling tij = t*j/Z 
is required. Such a scaling of the hopping amplitudes cannot be performed on 
the level of the Hamiltonian, but is possible in the effective action [139]. The 
bosonic DMFT ("B-DMFT") derived thereby treats normal and condensed 
bosons on equal footing and thus includes the effects caused by their dynamic 
coupling. The self-consistency equations of the B-DMFT are those of a single 
bosonic impurity coupled to two baths, one corresponding to bosons in the 
normal state and one to bosons in the condensate. The B-DMFT derived in 
the limit d or Z — > oo is again a comprehensive mean-field theory, i.e., it is 
valid for all input parameters and all temperatures. It not only reproduces all 
exactly solvable limits, such as the non- interacting (U = 0) and the atomic 
(tij = 0) limit, but also well-known approximation schemes for interacting 
bosons. For example, by neglecting all terms containing the hybridization 
function in the local action one obtains the mean-field theory of Fisher et 
al. [138]; for a detailed discussion see [139]. 

The solution of the self-consistent bosonic impurity problem defined by 
the B-DMFT equations requires new theoretical/computational methods. So 
far exact diagonalization [140, 141] and continuous-time QMC [142] were em- 
ployed. Thereby the bosonic Hubbard model was solved on the Bethe lattice 
for finite [140] and infinite [141] coordination number Z as well as on a simple- 
cubic lattice [142]. The phase diagram of correlated bosons on a simple-cubic 
lattice computed by the bosonic DMFT was found to agree with that obtained 
by numerically exact QMC to within 2% [142]. 

The B-DMFT is expected to be a valuable approximation scheme for the 
investigation of lattice bosons in situations where exact numerical computa- 
tions are difficult to perform or inefficient, as in the case of bosons with disor- 
der or many internal degrees of freedom, and for Bose- Fermi mixtures [143]. 
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1.8 DMFT for Nonequilibrium 

Recently the study of strongly correlated many-body systems out of equilib- 
rium has received much attention [144, 145]. This is motivated in particular 
by experimental progress in the investigation of ultra-cold atomic gases [133] 
and time-resolved pump-probe spectroscopy on strongly correlated materi- 
als [146, 147]. In general the real-time dynamics of correlated systems can be 
described by the extension of DMFT to nonequilibrium, provided that they 
are dominated by local temporal fluctuations and spatial correlations are not 
crucial. In nonequilibrium DMFT an effective impurity problem is formulated 
using the Keldysh formalism [148, 149], and as in equilibrium DMFT this 
mapping onto a single site becomes exact in the limit of infinite dimensions. 
For the Hubbard model in nonequilibrium the single-site DMFT action reads 

S = -i f dt H loc (t) - iJ2 I dt f dt' cl(t)A(t,t r )c a (t r ) , (1.49) 

where the Keldysh contour C runs from i m ; n to i max on the real time axis, 
back to i m i n , and finally to — i/3 along the imaginary time axis [150, 151]. 
The first term contains the local part of the Hamiltonian, e.g., H\ oc (t) = 
U(t)n^n^ — ^(n-j- + nj.) for a time-dependent interaction. The second term 
involves the hybridization function A(t,t') which couples the impurity to 
a time-dependent bath which must be determined self-consistently. Local 
contour-ordered correlation functions such as the Green function G(t, t') 
are obtained from the action (1.49) as expectation values (A(t)B(t r ) ■ ■ ■) = 
Tr[Tc exp(S)A(t)B(t') ■ ■ - ]/Z at the appropriate times, where Tq is the con- 
tour-ordering operator. For the Hubbard model this evaluation is the most 
demanding part of nonequilibrium DMFT and can so far be done with real- 
time quantum Monte Carlo methods [152, 151] for not too long times, and 
for sufficiently large U using a self-consistent perturbation expansion around 
the atomic limit [153]. For the Falicov-Kimball model, on the other hand, 
closed equations of motion govern the impurity Green function [19], which 
can be solved on the real time axis. The hybridization function A(t, t') in 
eq. (1.49) must be determined self-consistently by computing the local self- 
energy S(t,t') from the Dyson equation of the impurity model, calculating 
the momentum-dependent Green function of the lattice model from the lattice 
Dyson equation, integrating over momentum to obtain the local lattice Green 
function, and finally equating it with the impurity Green function. While this 
procedure is necessary for a Gaussian or other general DOS, for a semiellipti- 
cal DOS the self-energy can be eliminated and A expressed directly in terms 
of G [154]. 

Nonequilibrium DMFT was used to obtain the response of time-resolved 
photocmission [155-157] and optical spectroscopy [158] in correlated systems 
in terms of Green functions of the electronic system. The Falicov-Kimball 
and Hubbard models were studied in the presence of dc and ac electric fields 
[149,159-166,157], as well as for abrupt [167,158,156,152,151,153] or slow 
changes [168, 169] of the interaction parameter as a function of time. 

As an example, we consider a sudden change ( "quench" ) in the interation 
parameter of the Hubbard model from U = (i.e., with the non-interacting 
ground state as initial state) to finite values of U for times t > 0. For this case 
the DMFT equations were solved numerically for the paramagnetic phase and 
a semielliptic DOS [152, 151, 153]. In the following discussion the bandwidth 
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is equal to 4, and time is thus measured in units of 4fi/bandwidth, e.g., on the 
order of femto-seconds for a bandwidth on the order of eV. An interesting ques- 
tion is whether such an isolated system can thermalize due to the many-body 
interaction alone, i.e., whether at least some of its properties are the same as 
for an equilibrium system with the same energy [170]. Before the quench the 
momentum distribution is a step function, and the Fermi surface discontinuity 
An(t) remains nonzero for a finite time after the quench. For quenches to U 
< 2.5 the momentum discontinuity first reaches a so-called prethcrmalization 
plateau for t < 5 due to the vicinity of the integrable point at U = 0. This 
plateau in An(t) is given to good accuracy by 2Z — 1 [152], where Z is the 
Fermi-liquid quasiparticle weight in equilibrium at zero temperature and for 
interaction U. This value and also the transient behavior at short times is pre- 
cisely predicted by second-order unitary perturbation theory in U [171,172]. 
On the other hand, the double occupation essentially relaxes to its thermal 
value on this timescale, showing that the potential energy (and therefore also 
the kinetic energy) relax quicker than the occupation of individual states. For 
large U the behavior is different, showing strong so-called collapse-and-revival 
oscillations with approximate frequency 2tt/U. They stem from the vicinity 
of the atomic limit (i.e., zero hopping amplitude), for which the propagator 
e -tHt j g exac tly periodic with period 2-k/U [173]. For finite hopping (small 
compared to U) these oscillations are damped and decay on timescales of order 
ft/bandwidth. The oscillations of the double occupation are not centered at its 
thermal value, but rather at a different value that can be derived from strong- 
coupling perturbation theory [152]. The situation is thus similar to that at 
small coupling in the sense that the relaxation to the thermal state is delayed 
because the system is stuck in a metastable state close to an integrable point. 
However, both the weak-coupling prethcrmalization plateau in An(t) as well 
as the strong-coupling oscillations vanish in a narrow region of interaction pa- 
rameters U near 3.2 [152]. For quenches of U to approximately this value the 
system thermalizes rapidly: Both the momentum distribution and thus also 
the double occupation (due to energy conservation after quench) relax to their 
thermal values. In fact the retarded noncquilibrium Green function relaxes to 
the corresponding equilibrium function [151], so that all observables that can 
be obtained from it tend to the thermal value predicted by equilibrium sta- 
tistical mechanics. In other words, for quenches to interaction values in the 
vicinity of U « 3.2 the isolated system indeed thermalizes rapidly due to the 
many-body interactions. Many interesting questions remain open in this con- 
text, e.g., how thermalization depends on the parameters of the system. We 
refer to the reviews [144, 145] for further discussion. 

1.9 Summary and Outlook 

Due to the intensive international research over the last two decades the 
DMFT has quickly developed into a powerful method for the investigation 
of electronic systems with strong correlations. It provides a comprehensive, 
non-perturbative and thermodynamically consistent approximation scheme 
for the investigation of finite-dimensional systems (in particular for dimension 
d = 3), and is particularly useful for the study of problems where perturba- 
tive approaches are inapplicable. For this reason the DMFT has now become 
the standard mean-field theory for fermionic correlation problems, including 
cold atoms in optical lattices [174-176]. The study of models in nonequilib- 
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rium using an appropriate generalization of DMFT has become yet another 
fascinating new research area [149,151-153,156-169]. 

Until a few years ago research into correlated electron systems concen- 
trated on homogeneous bulk systems. DMFT investigations of systems with in- 
ternal or external inhomogcncities such as thin films and multi-layered nanos- 
tructures are still very new [177-182,175]. They are particularly important 
in view of the novel types of functionalities of such systems, which may have 
important applications in electronic devices. Here the DMFT and its non-local 
extensions [100-102] will certainly become very useful. 

In particular, the development of the ab initio band-structure calculation 
technique referred to as LDA+DMFT has proved to be a breakthrough in the 
investigation of electronically correlated materials. It provided already impor- 
tant insights into the spectral and magnetic properties of correlated electron 
materials [94,49,52,95-97]. Clearly, this approach has a great potential for 
further developments. Indeed, it is not hard to foresee that the LDA+DMFT 
framework will eventually develop into a comprehensive ab initio approach 
which is able to describe, and even predict, the properties of complex corre- 
lated materials. 
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